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PREFACE 


The work described herein was performed by the Control and Energy 
Conversion Division of the Jet Propulsion Laboratory. 
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DEFINITION OF SYMBOLS 

A surface area of critical nuclei particles, cm^ 

a inverse of the specific surface area of particles in fluidized 

beds, cm 

C concentration of gas phase, g raole/cm^ 

concentration of atomic hydrogen, g raole/cm3 
concentration of adsorbed atomic hydrogen, g mole/cm^ 
concentration of molecular hydrogen, g mole/cm3 
Ch2 concentration of adsorbed molecular hydrogen, g mole/cra^ 

Cg concentration of silicon vapor, g atom/cm3 
Cg concentration of adsorbed silicon, g abom/cra^ 

Cgjj concentration of silane in bulk vapor phase, g mole/cm^ 

Cgj^ concentration of silane at the substrate surface, g mole/cm3 

Cgg concentration of silane, g mole/cm3 

Cgj^ concentration of adsorbed silane, g mole/cra^ 

Cy concentration of vacant active sites, site/cm^ 

C[u concentration of total active sites, site/cm^ 

dp particle diameter, cm 

Dg surface diffusion coefficient, cm^/sec 

D diffusion coefficient, cra^/sec 

f a constant, 1 .554 X 10^ g1/2 cal''/2/cnj2 atm sec 

g* critical number of atoms (or molecules) in nucleus, g atoms 

(g) denotes gas phase 

G mass velocity, g/sec cra^ 

J rate of nucleation, no. of nuc^lei/cm^ sec 

Jp mass transfer coefficitirit . dim^snsionless 

k Boltzman constant 
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kg specific reaction rate constant for adsorption 

k(j specific reaction rate constant for desorption 

specific reaction rate constant for reverse desorption 
kh specific reaction rate constant for homogeneous pyrolysis 

kg specific reaction rate constant for surface reaction 

K. overall equilibrium constant for silicon CVD 
Kq mass transfer coefficient 

Kh equilibrium adsorption coefficient for atomic hydrogen 
Kh2 equilibrium adsorption coefficient for molecular hydrogen 

Kg equilibrium adsorption coefficient for silicon vapor 

Kgi equilibrium adsorption coefficient for silane 

m molecular mass, g/mole 

M molecular weight g/g mole 

n* concentration of critical nuclii, nuclii/cm^ 

n*g concentration of critical g-mers, g-raers/cm3 

ng concentration of g-mers, g-mers/cm^ 

nj^ concentration of gas, g mole/cm3 

Nh flux of hydrogen, g mole/sec cm^ 

Nr 

= ^si /(Ngi + Nh), flux ratio 
Nsi flux of silane, g raole/sec cm^ 

Ng average concentration of adsorbed silicon, g mole/cm^ 

Ng equilibrium concentration of adsorbed silicon, g mole/cm^ 

p partial pressure of gas, atm 

p' equilibrium vapor pressure of gas, atm 

Pgi partial pressure of silane, atm 

r* radius of the critical nucleus, cm 

rg rate of adsorption, g mole/sec cm2 


V 
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Pjj rate of desorption, g mole/sec cm2 

rh rate of homogeneous pyrolysis, g mole/sec cm^ 

r^p rate of heterogeneous pyrolysis, g mole/sec cm3 

rjn rate of mass transfer, g mole/sec cm^ 

rQ overall rate of silicon CVD, g mole/sec cm3 

rs rate of surface reaction, g mole/sec cm^ 

R gas constant 

Rep particle Reynolds number, dpG/u, dimensionless 
(s) denotes solid phase 

S supersaturation ratio, p/p' 

Sh Sherwood number, dimensionless 

T temperature, K 

Tg substrate temperature, K 

V reactor volume 

X root mean square distance for diffusion before reevaporation 

(or desorption), cm 

Z, film thickness for gas phase diffusion, cm 

z' Zeldovich correction factor 

ag condensation coefficient 

Af vapor phase free energy change, cal/g 

AF* critical free energy change for nucleation, cal/nucleus 

AFb bulk free energy change for nucleation, cal/nucleus 

AFg surface free energy change for nucleation, cal/nucleus 

AFg free energy change for nucleation as a function of number of 
molecules in the nucleus, cal/g mole 

AF* Gibbs free energy change of activation for surface diffusion 
cal/g mole 

AF* Gibbs free energy change of activation for desorption, cal/g mole 
des 
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AF °2 standai"d free energy of hydrogen molecules, cal/g mole 

AF° standard free energy of silicon, cal/g atom 

s 

AF°i standard free energy of siiane, cal/g mole 

AF° standard free energy change for reaction , cal/g mole 

AH° standard heat of formation, cal/g mole 
AH° standard heat of reaction, cal/g mole 
e void fraction 

6t eddy diffusivity coefficient, cm^/sec 

average eddy diffusivity coefficient, cm^/sec 
6 contact angle , deg 

X mean spacing between ledges, cm 

viscosity coefficient, g/cm sec 
V surface vibration frequency, sec"^ 

p density, g/cm3 

cr specific surface free energy, erg/cm^ 

‘■'■g_v surface tension between substrate-vapor, erg/cm2 

surface tension between condensate-substrate, erg/cm^ 
<^_v surface tension between condensate-vapor, erg/cm^ 

cj) bulk flow correction factor 

X functional relationship for contact-angle 

molecular volume 

o) nuclei growth frequency, molecules/nucleus sec 
CO frequency factor for activated species 
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ABSTRACT 


The four basic elements in the chemical vapor deposition (CVD) 
of silicon from silane are mass transport of silane, pyrolysis of silane, 
nucleation of silicon and silicon crystal growth. These four elements 
are analytically treated from a kinetic standpoint. Rate expressions 
that describe the various conceivable steps involved in the chemical 
vapor deposition of silicon are derived from elementary principles. 
Applications of the rate expressions for ( 1 ) modeling and the simulation 
of the silicon CVD process and (2) the analysis of experimental data 
on silicon CVD are discussed. The lack of an experimentally established 
mechanism of the silicon CVD process and established values for var- 
ious constants involved in the rate expressions is the major impediment 
to the modeling of the CVD process. Experimental data are needed to 
determine the equilibrium adsorption coefficients for silane, hydrogen 
and silicon vapor and the " .tivation energies and frequency factors 
for the various rate proce ^es involved in the silicon CVD. 
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SECTION I 
INTRODUCTION 


The Low-Cost Silicon Solar Array (LSSA) project sponsored at JPL 
by ERDA has the responsibility for developing the technology and the 
production capability for economically practical solar photovoltaic 
electricity generating plants. The objective is the achievement of a 
generating capacity of 500 megawatts peak at a cost of $500 per kW by 
1985. The program of LSSA is divided into four technology tasks, each 
of which has the responsibility for advancing a particular technical 
area so as to achieve the production and cost goals, 'i’he first of 
these tasks is the Silicon Material Task, for which the 1986 objective 
is to reduce the price of Si material, suitable for- the fabrication 
of solar cells having satisfactory performance, from the present price 
of about $65/kg (for semiconductor grade Si) to less than $ 10/kg. 

To achieve the price objective of the Si Material Task, new areas 
of chemical processing and solar cell technology must be studied and 
developed, since no modification of the basic processes presently used 
commercially will lead to the attainment of the $ 10/kg price. The 
inability of the present commercial process to meet the task objective 
is revealed from a consideration of the features of this process. 

The first step is the reaction ,">f metallurgical grade Si, having a 
purity of about 98.5$, with HCJ. in a fluidized bed reactor to yield 
the intermediate SiHCl3 in a mixture of chlorinated silanes. After 
chemical treatment and fractional distillation to obtain very pure 
SiHCl3, semiconductor grade Si is produced by chemical vapor deposition 
resulting from the H2 reduction of SiHCl3. This process to produce 
extremely pure Si is costly and energy-intensive, the final deposition 
consuming. an estimated 385 kWh/kg Si product. 

Several different chemical processes are being investigated under 
JPL contracts in the task program. A contract with Union Carbide involves 
the hydrogenation of SiClij to SiHCl3 in a fluidized bed as well as 
a process for the deposition of Si from SiH4 in another fluidized bed. 

To support the developments in silicon production, the task program 
includes an in-house JPL subtask for supporting studies of silicon pro- 
duction as well as a subtask for consultation in this area by Professors 
0, Levenspiel and T. Fitzgerald of Oregon State University. The in-house 
program consists of efforts in modeling of the deposition of Si in a 
fluidized bed reactor, experimental studies of SiHi^ pyrolysis and fine 
particle fluidization, and thermodynamic and chemical engineering analysis 
of the reactions in the fluidized process for silicon production. 

The initial phase of the effort in analysis of silicon chemical vapor 
deposition is described in this report. 

The finite amount of available fossil fuels and the continually 
increasing energy requirements of mankind dictate the need for the 
development of alternate energy sources. Solar energy is one energy 
source that is abundantly available and can be harnessed by silicon 
solar cells to generate electric power. The production of polycrystalline 
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silicon by the pyrolysis of silane is a promising method for reducing 
the present-day high cost of silicon material suitable for silicon 
cells. Although the pyrolysis of silane has been extensively studied 
in connection with epitaxial growth of single-crystal silicon films, 
the mechanisms of .silt 'Ion deposition in fluidized bed and free space 
reactors are little understood. The deposition process consists of 
several physical and chemical steps. The controlling step will depend 
on the reactor temperature and silane concentration. It is the purpose 
of this paper to establish the thermodynamic feasibility of the overall 
pyrolysis reaction and to develop rate expressions for the various 
individual steps of the deposition process. Rate expressions for the 
heterogeneous reaction kinetics of silane pyrolysis will be developed 
in the fashion of Hougen and Watson's treatment of Langmuir-Hinschelwood 
theory. Excellent accounts of this type of analysis can be found in 
Hougen and Watson (Reference 1-1) and Smith (Reference 1-2). Hirth 
and Pound (Reference 1-3) and Powell et al. (Reference 1-4) are the 
references used most in developing the kinetic expressions for nucleation 
and growth of silicon crystals. 


Experimental data on the pyrolysis of silane have been reported 
by many investigators (R ‘ences 1-5 through 1-9) usually in the form 
of Arrhenius equations, i.^uiietic parameters (activation energy and 
frequency factor) are studied for the calculation of the specific reaction 
rata constants for the conversion of silane or growth rate of silicon 
crystals. Different experimental conditions (temperature, pressure, 
concentrations, flow rates, and reactor geometry) have been used, and 
hence different overall kinetic constants for the pyrolysis of silane 
have been reported. These differences are generally attributed to 
the differences in the mechanisms of the CVD process under different 
experimental conditions, although the exact mechanism has never been 
proved. 


Experimental evidence leads to some general conclusions about 
silane pyrolysis; the reaction is predominantly homogeneous at low 
temperatures (300-500° C) and is predominantly heterogeneous at higher 
temperatures (500-1200°C) . The heterogeneous reaction is a diffusion 
limited process above 1000°C. The crysta.1 structure of deposited silicon 
is a strong function of substrate temperature. If the reaction rate 
constants and expressions are known for all the individual steps in 
all possible mechanism sequences, it will be invaluable in the calculation 
of the overall conversion or growth rate for specified experimental 
conditions. 


Development of rigorous rate expressions is the essential first 
step in this direction. In the present analysis, no assumptions are 
made as to the mechanism of the silicon CVD process, and rate expressions 
are derived for every conceivable elementary step (within the limits 
of reason). It is hoped that the present analysis will provide analytical 
expressions for modeling of the silicon CVD process and provide a sound 
theoretical basis for future experimental research on silicon CVD. 
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SECTION II 

THERMODYNAMICS SILANE PYROLYSIS 


There are two powerful tools in thermodynamics for the analysis 
of any chemical reaction, viz: the heat of reaction and the free energy 
change for the reaction. These two quantitative measures determine 
the feasibility of the reaction, the maximum extent of conversion and 
energy requirements for the process. The heat of reaction and free 
energy of reaction can be calculated from the heats of formation and 
standard free energies of the reactants and products by the following 
relations : 


heat of reaction (AH°) = heat of formation ~ heat of formation of 

^ of products (AH°) reactants (AH° ) 


standard free 
energy of 

reaction (AF° ) 
R 


standard free 
= energy of 

silicon (AF° ) 
s 


standard free 
+ energy of 

hydrogen (2AF° ) 
H2 


standard free 

energy of 

silane (AF° ) 
si 


The reaction equilibrium constant can be obtained from the expression 


(AF°)r = - RT In k 


The overall reaction for silane pyrolysis can be written as 


SiH^Cg) ► Si(s) + 2H2(g) 

The standard free energy change, equilibrium constant and heat 
of reaction have been calculated for the temperature range 100 - 1600 K 
and are tabulated in Table 2-1 . As can be observed from the data in 
Table 2-1 , the pyrolysis of silane is exothermic and thermodynamically 
feasible in the temperature range 100 - 1600 K. The high values of 
the equilibrium constant, especially at higher temperatures, indicate 
almost 100$ conversion, hence the irreversible nature of the reaction. 
Thus the thermodynamic analysis has shown that pyrolysis of silane 
is an irreversible and exothermic reaction. 
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Table 2-1. Standard Free Energy Change, Equilibrium Constant, 
and Heat of Reaction 

Reaction: SiH^^Cg) »-Si(s) + 2H2(g) 

units of free energy change and heat of reaction: cal/mole 


Temp, K AFg^ AF| ^^2h2 (-^r) In K -AH| 


100 

200 

298.15 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 


3164.4 

7433.1 

12056.4 

12146.9 

17206.6 

22574.5 

28225.6 

34139.4 

40298.4 

46686.8 

53290.4 

60096.1 

67092.1 

74267.6 

81612.7 

89118.2 

96776.0 


27.9 

212.2 

571.8 

580.2 

1108.3 

1769.6 

2544.0 

3415.9 

4374.0 

5407.5 

6510.8 

7677.7 

8903.3 

10183.6 

11515.4 

12893.5 

14318.9 


3359.6 

8685 

14561.6 

14678.0 

21138.4 

27949.4 

35042.2 

42371 .0 

49900.6 

57606.2 

65469.4 

73482.8 

81626 

89891 .8 

98280.4 

106770.8 

115362.8 


223.1 

1464.1 

3077.0 
3111.3 

5040.1 

7144.5 

9360.6 

11647.5 

13976.2 

16326.9 

18689.8 

21064.4 

23437.6 

25807.8 

28183.1 

30546.1 

32905.7 


+1.1277 

+3.6841 

+5.1033 

5.2189 

6.3407 

+ 7.1905 

7.8507 

8.3732 

8.7913 

9.1289 

+9.4050 

9.6364 

9.8285 

9.9900 

10.1302 

10.2476 

10.3492 


10097.7 

9311.1 

8634.9 

7800.0 

7016.6 

6394.6 

5917.9 

5569.7 

5330.1 

5177.2 

5091.9 

5057.7 

5060.5 

5088.0 

5132.2 

5184.5 

5240.7 


Heat of formation and free energy data are from McBride, Heimel, 
Ehlers and Gordon, Thermodynamic Properties to 6000° K for 210 
Substances Involving the First 18 Elements . NASA SP-3001, 1963. 
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SECTION III 

KINETICS OF SILICON CVD 


The overall chemical vapor deposition process consists of several 
physical and chemical steps. All the possible steps are illustrated 
in Figure 3-1 • The exact path of the deposition mechanism depends 
on the temperatures of the vapor phase and substrate and silane 
concentrations . 

The conceivable elementary steps* involved in silicon CVD are: 

(R-l) homogeneous pyrolysis of silane 

homogeneous nucleation of silicon vapor 
(R 3 ) mass transport of silane gas to the substrate surface 
(R 4 ) mass transport of hydrogen away from the substrate surface 
(R 5 ) heterogeneous pyrolysis of silane adsorbed directly on nuclei 
( 85 ) adsorption of silane on substrate surface 
(R 7 ) adsorption of silicon vapor on substrate surface 
(Rg) desorption of hydrogen from substrate surface 



Figure 3-1 . Schematic of Silane Pyrolysis 


These are the steps that are assumed to prevail , as the exact mechanisms 
of silicon CVD are not known. 
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(Rg) heterogeneous pyrolysis of silane adsorbed on substrate 
surface 

(R-IO) heterogeneous nucleation of adsorbed silicon 

These elementary steps can be conveniently grouped into six cate- 
gories: (1) mass transport of silane and hydrogen, (2) homogeneous 

pyrolysis of silane, (3) homogeneous nucleation of silicon, (4) heter- 
ogeneous pyrolysis of silane, (5) heterogeneous nicleation of silicon, 
and (6) silicon crystal growtn. The actual mechanism of silicon CVD 
will consist of three or more of these six categories. Kinetic rate 
expressions for each of these steps will now be developed. 
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SECTION IV 
MASS TRANSPORT 


The diffusion of the reactant (silane) to the substrate surface 
(step Rg) and the diffusion of the product (hydrogen) away from the 
substrate surface (step R 4 ) are the important mass transport processes 
involved in the chemical vapor deposition of silicon from silane. 

The rates of R^ or Rij are not independent of each other and are dependent 
on the stoichiometry of the pyrolysis reaction. In the absence of 
equilibrium limitations, diffusion can be the limiting process in the 
CVD uf silicon at high temperatures. In fact, many inve*stigators of 
silicon epitaxy reported the process to be diffusion-limited. In such 
instances the reaction kinetics will be extremely fast and the overall 
rate of CVD can be obtained from a knowledge of the mass transport 
rates. The rate expressions for mass transport in flow reactors can 
be written by using the Stephen-Maxwell relation as; 


N 


si 


8Cgs 

= _(D + € .) + 

dZ 



(N 


SI 


Nh) 


( 1 ) 


where Ngj_ is the flux of "dlane, D is the binary diffusion coefficient, 
is the eddy diffusivity coefficient, Cgg is the concentration of 
silane, C is the concentration of total gas phase, N^ is the flux of 
hydrogen and Z is the diffusion path length. Molecular diffusion is 
accounted for by D, and is an analogous eddy diffusivity to account 
for mass transport by turbulent convection. Eddy diffusivity is a 
function of distance from the surface, and this relationship is not 
generally known. Thus an average constant eddy diffusivity is used 
for integration of Eq. (1) from bulk silane concentration (Cgb) to 
the silane concentration at the surface (Cg^). The result of integration 
is: 



The value of Z, the diffusion path is not generally known and, 
hence the term 

r-p) 
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will be replaced by one variable, K^, called a mass transfer coefficient. 
Equation (2) can be simplied by defining Np and cj) : 


Nr = 


’SI 


Nsi + Nh 


<P = 


f ^ss\ 

- Tr 

Nr 


(3) 


(4) 


with LM denoting logarithmic mean. For silane pyrolysis, 2Ng;j^ = 
by reaction stoichiometry and hence Nr = -1 and <|> becomes: 


*^si ^sb 



Equation (2) can now be reduced to a more familiar form of mass transfer 
equations, where r,j] is the rate of mass transfer per unit reacter volume 
and a is the inverse of the specific surface area of particles in the 
reactor. 

Kc 

^'m ~ “T ^ (*^sb - *^si^ 

9 

But, as shown in Eq. (5), <i’ is a function of Cgj^ and Cg^ and cannot be 
included in the mass transfer coefficient. Substituting the value of 
4) from Eq. (5) in Eq. (6), the rate of mass transfer in silane pyrolysis 
reactors can be expressed as 


r,n = k(, a C In 


1 + 


Csb 


1 + 


^31 

T j 
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For fluidized bed reactor systems the value of a' is expressed as 


, 6(1 - O 

a 

dp 


f 


where e is the void fraction in the fluidized bed and dp is the diameter 
of the particles in the reactor. The value of Kg can be obtained from 
empirical correlations. However, as ^leankoplis (Reference 4-1) points 
out, care should be exercised in using these correlations, as some 
of the experiments were not conducted with <)) = 1 . In such instances, 
the value of for the experimental conditions used in developing the 
correlations should be calculated and the true value of Kq obtained 
therefrom. 

If the mass transport is the rate-limiting step in the silicon 
CVD process, i.e., the mass transport is the slowest step, the concen- 
tration of silane at the substrate can be assumed to be zero (Cg^ = 0). 
Then, the rate of silicon CVD process can be determined from Eq. (7). 
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SECTION V 

HETEROGENEOUS PYROLYSIS 


The traditional approach to the analysis of heterogeneous kinetics, 
utilizing Langmuir's theory of adsorption and Hinschelwood 's theory of 
surface reactions, was developed in detail by Hougen and Watson (Refer- 
ence 1-1). The major assumptions made in this approach are (1) the 
reaction surface is energetically uniform, i.e., all the active centers 
have the same activity for adsorption; .(2) the amount adsorbed has no 
effect on the rate of adsorption per site; i.e., there is no interaction 
between adsorbed molecules; (3) all the adsorption occurs by the same 
mechanism and each adsorbed complex has the same structure; (i|) the 
extent of adsorption is less than one complete raonomolecular layer on 
the surface; (5) the total number of active sites is constant; and 
(6) one of the steps in the process is usually controlling and the 
rest occur at equilibrium. 

The three steps involved in heterogeneous decomposition of silane 
are adsorption of silane on the silicon substrate (R6), surface reaction 
of silane (Rg), and desorption of hydrogen (Rs)- Rate expressions for 
the three steps are derived by the Hougen-Watson procedure and included 
in Table 5-1. Rate equations for the surface reaction and hydrogen 
desorption are different for the different mechanisms of hydrogen adsorp- 
tion (molecular, atomic or none). Stoichiometric equations and equilib- 
rium concentrations of adsorbed species are also included in Table 5-1. 
Overall rate expressions for heterogeneous pyrolysis of silane can 
easily be derived if it is assumed that one of the steps is rate con- 
trolling and the rest occur at equilibrium. One useful relation in 
deriving the overall rates is Eq. (9), in which ’0^’ represents the con- 
centration of tota^ active site£, 'C represents the concentration of 
vacant sites, and Cg^, and Cg represent concentrations of adsorbed 
silane, hydrogen and silicon: 


^m " ^v * ^si + + Cg 


(9) 


The overall rates (r^p) of heterogeneous pyrolysis of silane are 
derived using Eq. (9) and rate equations and equilibrium concentrations 
given in Table 5-1. Equations (10) through (17) are the rate equations 
derived with one step assumed as the controlling step. 

Silane adsorption controlling, with molecular hydrogen adsorption: 


f’hp = 


^a^m 
1 + Kv,oC 


h2^h2 


^31 


( 10 ) 
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Table 5-1. Rate Expressions for Heterogeneous Pyrolysis of Silane 




Hydrogen Adsorption Mode 


Controlling Step 

Molecular 

Atomic 

None 

Adeorption of eilane 

SiHp(g) + X -SIH4 • X 

/ CgA 

•'a ■ •'af’^si^v “ 

SiH4(g) + X SiH4 • X 

/ Csi\ 

■'a ■ *^a(^si^v “ j 

SiH4(g) + X SiH4 • X 

/ Cgi\ 

'’a ■ *^al*^3i*^v ~ ) 


(Caileq. = •^si*'3i^v 

(Cgi)eq. = f^si^si^'v 

(Csi)eq. = 


Surface reaction 


SiHis . X + 2X Si • X + 2 H 2 X 


/ Cv\2 


-3 = >^3Csil - 


(^3i^eq. = ^ 


SiHij . X + 4X — -Si • X + 4H X SiH4 • X -Si • X + 2 H 2 


/ Cv\4 


*'s - '^3^3il ~ 


( C3i ) eq . = 0 




( Cgi ) eq . = 0 


De3orption of hydrogen r^ = -k(j(CvCh2 


H2 • X -H2(g) + X 

Ch2 

^^h2^eq. ~ *^h2^v*^h2 


2H • X H2(g) + 2X 

-kd/ 2 Ch 

= srr'' ■ 


=■2 


^^h^eq. ■ *^h*^h2^v^ 
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Surface reaction controlling, with molecular hydrogen adsorption; 


^s^si^m 


•"hp “ 


(1 + Ksj_Cgj_ + Kh2Ch2)^ 


■'SI 


( 11 ) 


Hydrogen desorption controlling, with molecular hydrogen adsorption: 


*'hp = 




( 12 ) 


^h2 


Silane adsorption controlling, with atomic hydrogen adsorption: 


^hp = 


^a^ra 


1 + \/Kh^ 


'SX 


( 13 ) 


Surface reaction controlling, with atomic hydrogen adsorption: 


*'hp = 


^s^si^m 


(1 + KgiCgj_ + x/KhCh2>^ 


'SI 


( 14 ) 


Hydrogen desorption controlling, with atomic hydrogen adsorption: 


^d^m 


•’hp 


( 15 ) 


Silane adsorption controlling, with no hydrogen adsorption: 


*’hp " ^a*^ra^si 


( 16 ) 


Surface reaction controlling, with no hydrogen adsorption: 


•’hp " 


^s^si^m 

1 + KsiCgj^ 


'SI 


( 17 ) 
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The overall rate of heterogeneous decomposition of silane can 
be described by one of the equations (10) - (17), depending on the 
controlling step. If adsorption, surface reaction and desorption are 
all important, then overall rate at steady state conditions is given 
by the relation: 


I'O = ‘'a = »*s = ^d 


The concentrations of adsorbed silane, silicon and hydrogen can 
be evaluated from the three independent equations, implicit in the 
above relation. Thus the overall rate of heterogeneous pyrolysis of 
sij.ane can be expressed in terms of gas phase concentrations and the 
equilibrium adsorption coefficients. The equilibrium adsorption coef- 
ficients can be evaluated from independent adsorption studies. Since 
the equilibrium adsorption constants and specific rate coefficients 
are all functions of temperature, a large amount of kinetic data at 
different temperatures will be necessary to evaluate all the equilibrium 
and rate constants from kinetic data alone. An example of this type of 
analysis is given by Hougen and Watson (Reference 1-1) for the hydrogena- 
tion of Codimer. As had been pointed out by many investigators, it may 
be extremely difficult, if not impossible, to determine the mechanism 
and identify the rate-limiting step from kinetic data alone. Experimental 
kinetic data may fit equally well to more than one overall rate expression, 
if all the equilibrium and rate constants are derived from the kinetic 
data. Thus it is very important to obtain the equilibrium adsorption 
constant from independent adsorption studies, if the rate equations 
are to be used with confidence in any modeling effort. 

The kinetic parameters studies of activation energy and frequency 
for the overall heterogeneous pyrolysis of silane reported in the literature 
vary greatly. Joyce and Bradley (Reference 1-9) reported an activation 
energy of 37 kcal/mole for the epitaxial CVD of silicon on silicon sub- 
trates with substrate temperatures in the range of 920 - 1260°C. They 
indicated that mass transfer is the rate-limiting step in the temperature 
range of 1100 - 1260°C and that either chemisorption, surface reection 
or desorption is the limiting step in the temperature range of 950 

- 1100°C. Their experiments were carried out at low silane pressures 
(2 mm Hg) to suppress homogeneous decomposition. They speculated that 
surface reaction or desorption of hydrogen is the rate-limiting step 
in the low-temperature (950 - 1100°C) range. Henderson and Helm 
(Reference 1-8) reported an activation energy of 20 kcal/mole for the 
growth of homoexitaxial growth of silicon on silicon substrates using 
silane pyrolysis at temperatures of 823 - 983°C and pressures of 0.15 

- 0.02 torr. From an analysis of their data they surmised that surface 
chemical reaction is the controlling step, silane is weakly adsorbed 
and hydrogen is not adsorbed. Farrow (Reference 1-7) reported an acti- 
vation energy of 17 kcal/mole for CVD of silicon on silicon substrates 
and that his data agreed well with the reaction mechanism proposed by 
Henderson and Helm (Reference 1-8). 

Everstiyn and Put (Reference 1-6) conducted the CVD of silicon 
on silicon substrates from the pyrolysis of silane at atmospheric pressure 
and a temperature range of 590 - 1050°C. They reported an activation 
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energy of 37 kcal/mole. The differences in the values of activation 
energy reported are generally attributed to impurities on the substrates 
and different surface orientations. The rate limiting step has not been 
established clearly in the literature. Experimental data at atmospheric 
pressure are also very limited. It therefore seems imperative to carry 
out an experimental program to ( 1 ) obtain kinetic data on the CVD of 
silicon from silane pyrolysis at atmospheric pressure and (2) establish 
the rate-limiting step for various experimental conditions. 
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SECTION VI 

HOMOGENEOUS PYROLYSIS 


Homogeneous decomposition of silane (step Ri in Figure 2-1) is 
a simple first-order vapor phase pyrolysis of silane. The stoichiometric 
equation can be written as: 


SiH4(g) ►Si(g) + 2H2(g) 

The rate of reaction r^ can be expressed as 




*^h*^si 


( 18 ) 


Joyce and Bradley (Reference 1-9) reported an activation energy of 
51 kcal/mol. Everstiyn (Reference 1-5) reported a value of 30 kcal/mol 
after conducting extensive research on vapor phase pyrolysis of silane 
in the temperature range of 660 - 755“C. Eversteyn determined also 
the critical concentrations of silane necessary for homogeneous reaction 
in the temperature range 780 - 1136°C. This critical concentration 
decreased from a value of 0.248H at 780°C to a value of 0.002il at 1136°C. 
The logarithm of the critical concentration varies linearly with the 
inverse of the temperature with an activation energy of 38 kcal/mol. 

Joyce and Bradley (Reference 1-9) reported that homogeneous silane 
pyrolysis can be totally suppressed if the silane pressure is kept 
below 2 ram Hg. 
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SECTION VII 

HOMOGENEOUS NUCLEATION 


In the absence of all foreign materials and surface, vapor 
condensation is blocked by an activation free energy barrier which 
arises from surface free energy increases resulting from appearance 
of embryos of the more condensed phase. This can be seen more clearly 
through free energy considerations. In homogeneous nucleation processes 
(step R 2 in Figure 3-1) it is necessary to go to higher super saturations 
to get over the free energy barrier. The excess free energy of a vapor 
at a supersaturated pressure p over that of a liquid at the same temper- 
ature T is given by the relation 


Af = RT In (p/p') 


where Af is the excess free energy per gram, R is the gas constant 
per gram of vapor, and p is the saturation vapor pressure at temperature 
T. This Af will be the decrease in the system free energy as an embryonic 
droplet is formed. At the same time there will be a free energy increase 
associated with the embryo's surface free energy, which can be written 
as (4/3)Trr3o-, where r is the radius of the embryo and cr is the specific 
surface free energy (ergs/cm^). The phase transition can occur only 
through drop-wise condensation (as all foreign materials are absent). 

Thus the free energy decrease by phase transition is dependent on the 
volume of the embryo and can be written as (4/3)'n'r3 PRT In S, where 
P is the density of the condensed phase and S is the supersaturation 
ratio (p/p ). Thus the net free energy change for embryo formation 
can be shown to be 


AF = 47rr2a- - (4/3)n-r3pRT In S 


(19) 


The free energy change has to be negative for successful vapor condensation. 
At supersaturation ratios of 1 and <1 AF is always positive and 
phase transfer does not occur. However, as S increases, AF goes through 
a maximum. This is the maximum free energy barrier (AF*) to be surmounted 
for vapor condensation (Figure 7-1). 

Invoking the basic principle that any process in which AF decreases 
is a spontaneous process, we can see that if, by any .means, the embryo 
can reach the maximum AF, the slightest additional condensational growth 
of the embryo beyond the size associated with the maximum AF should 
lead to irreversible growth of the embryo from supersaturated vapor. 

Thus the radius of the embryo r for which AF attains its maximum will 
be a critical embryo radius for whatever S value present. This critical 
radius can be obtained by differentiating Eq. (19) with respect to 
r at constant T and S and setting the 9(AF)/9r equal to zero and then 
solving for r* (critical radius). 
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Figure 7-1 . Free Energy of Formation of Jritical Nuclei 


r“ = 20-/PRT In S 


( 20 ) 


This r* is the radius of the embryo that is essential for vapor conden- 
sation at temperature T and supersaturation S. Thus r* is the smallest 
diameter of the condensed phase. The g*, the number of molecules or 
atoms in an embryo of critical radius r*, is given by the relation 
((4/3)'n’r3p = gm) where m is the mass of 1 molecule, 


S* 


32ir cr 3 

3mp2R3x3(in s)3 


( 21 ) 


So far, consideration of the thermodynamic aspects of nucleation phenomena 
has created a paradox. The supersaturated vapor is in a raetastable 
state because its specific free energy exceeds, by the finite amount, 

RT In S, that of condensed phase at the same temperature. However, 
it cannot jump to the preferred lower free energy state because that 
jump can only proceed via formation of tiny drops, beginning with clusters 
of two, three, four molecules, etc. Th" growth of such embryos of the 
new phase carries the system up to states where the total free energy of 
the system has risen more by creation of new surface free energy of the 
embryonic drops than it has fallen by virtue of the bulk free energy 
change accompanying the phase transition that produced these drops. 

That is, phase transition is always, in the homogeneous case, blocked 
by an activation energy barrier. 

This barrier can be overcome only by raising S high enough, which 
will lower the AF* and the r* that must be exceeded in order that the 
system will spontaneously shift bo the preferred condensed phase. 

Even at subsaturations (S<1), there exists, as a steady-state population, 
a distribution of embryos satisfying a Boltzmann-type relation: 

ng = n^e-^S/i^T (22) 
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in which ng is the number of clusters of g molecules and AFg is in 
terms of ’g’; k is the Boltzmann constant. If rare fluctuation phenomena 
can, even at S<1, support such a population of molecular clusters, 
then when S rises above unity, surely a general increase in that popula- 
tion is expected at all g. Thus, for some large S (Sw5,6) it can 
be expected to reach a point where fluctuations can send some embryos 
over the top of the energy barrier whereupon the supercritical embryos 
can subsequently grow without limit until the vapor pressure is pulled 
from p down to p and the system is thermodynamically stabilized. 
Substituting Eq. (20) for r* in Eq. (19), the free energy barrier to 
be surmounted for the formation of critical embryos of size r* is 
obtained: 


16tt ct3 

AF» (23) 

3(pRT In S)2 


Hirth and Pound (Reference 1-3) enumerated the classical nucleation 
theory of Volmer and Weber, in which the rate of nucleation is formulated 
in terms of equilibrium population critical embryos ng, surface area 
of embryos A and an impingement frequency factor gj. Volmer and Weber's 
equation for homogeneous nucleation rate is 


J = Aojn» (2^0 

g 

where A = AiTr*2, u = aQp/(2iTmkT) and ng# is given by Eq. (22). 
Substituting these expressions in Eq. (24;, we obtain 


J 


ac(4Trr*2) 


P 

(2iTmkT)1/2 


me-AF*/kT 


(25) 


where is the condensation coefficient for molecules impinging on 
the critical embryos. 


Hirth and Pound (Reference 1-3) also described a refinement of 
Volmer and Weber's theory by Becker and Doering and Zeldovich. In the 
refined treatment, instead of using equilibrium concentrations the 
actual steady-state concentrations of critical nuclei were calculated 
by considering the changes in cluster size by the addition and removal 
of a molecule. This elegant approach resulted in the incorporation 
of an additional term Z’, the so called nonequilibrium factor, into 
Eq. (24). The refined expression for homogeneous nucleation rate then 
becomes , 

J = UqZ ' (4Trr*2 ) 


(i;TrmkT)^/2 


n-] exp(^F*/kT) 


(26) 
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where 


, m (PRT In S)^ 

Z 

P 8 ■n'tr(crkT)^^^ 


Substituting Eq. (20) for r* and noting that p = nk:‘, Eq. (26) can 
be simplified to yield 


J 



exp(-AF*/kT)n^ 


(27) 


McDonald (Reference 7-1), in a similar treatment, considering the 
rate of growth and decay of critical embryos, arrived at a homogeneous 
nucleation rate expression identical to Eq. (27). The Inherent fallacy 
in Eq. (27) arises from assigning macroscopic thermodynamic properties 
to microscopic embryos. Many corrections have been proposed in the lit- 
erature to correct this problem. For this reason, for critical embryos 
in the size range 1 < g* <50 a different statistical mechanical model should 
be used in place of Eq. (27). Lothe and Pound (Reference 7-2) suggested 
another correction factor to Eq. (27) to take into account the contribu- 
tions of molecular entropy of mixing of small clusters with monomer, 
rotational and vibrational free energies of initial embryos to the free 
energy of formation of critical embryos. However, all these corrections 
can be taken care of by assigning a value sufficiently less than unity 
to the condensation coefficient a,,. Equation (27) is strictly applicable 
to nucleation of droplets from unary vapors. It had been reported in 
the literature that the same equation can be applied to nucleation of 
crystallites from both unary and binary vapor with proper modifications 
to a,, or effective supersaturation ratio. 
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SECTION VIII 

HETEROGENEOUS NUCLEATION 


Heterogeneous nucleation (step R-|o in Figure 2-1) is the process 
by which the adsorbed silicon atoms on the substrate are incorporated 
in the silicon deposit. The thermodynamics of heterogeneous nucleation 
is vary similar to that of homogeneous nucleation. The added parameter 
is the equilibrium contact angle of the nuclei . A schematic model 
of nucleation of a spherical cap-shaped nucleus on a solid substrate 
is shown in Figure 8-1 . 


The free energy of formation of a spherical cap-shaped nucleus 
consists of two terms; the bulk free energy change AFj^ and surface 
free energy change AFg. 


AFb = -Trr3A'(e)Af 


( 28 ) 


Af = -pRT In S 


(29) 


X(0) = 


(2 + cos 0) (1 - cos 0)2 

- 


(30) 


where A'(G) is a geometric factor relating the volume of a spherical 
cap to that of a sphere of radius r, and 0 is the equilibrium contact 
angle determined by the surface free energy equilibrium. 


°"s-v - *^c-s + *^c-v cos 0 (31) 



Figure 8-1. Spherical Cap-Shaped Nucleus 
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and where cTg_y, t^'c-s o'g.y are the specific interfacial free energies 
of the substrate vapor, condensate-substrate, and condensate-vapor 
interfacing, respectively. The surface free energy term is given by 

AFg = 4Trr2x(e)o- (32) 


Therefore, the total free energy change associated with the formation 
of nuclei is 


AF = k-nr^X (0)tr 



X(0)pRT in S 


(33) 


The variation of the free energy change with r is very similar to that 
for the homogeneous nucleation, as shown in Figure 7-1, and passes 
through a maximum corresponding to the critical nuclei radius r*. 

The value of this r* can be found by maximizing AF with respect to 
r, and r* is found to be: 


and 


2cr 

V* = 

pRT In S 


AF» 


16 y„3x(0) 

3 (pRT In S)2 


(34) 


(35) 


By comparing Eqs. (34) and (35) with Eqs. (20) and (23) > it can be 
readily observed that the radius of critical nucleus is the same in 
both homogeneous and heterogeneous cases, but the critical free energy 
change in the heterogeneous case is dependent on the contact angle 0 . 

If the contact angle 0 is zero (X(0) = 0, complete wetting of the sub- 
strate by the condensate) , AF* is zero and nucleation is most rapid. 

If 0 is equal to 180 deg, corresponding to no wetting of wall at all, 

X(.B) = 1 and Eq. (35) reduces to Eq. (23), the expression for homogeneous 
nucleation. This is the limiting case in which the substrate has negli- 
gible potency as a nucleation catalyst. 

Heterogeneous nucleation kinetics are formulated much the same 
way as homogensous nucleation kinetics and the rate of nucleation is 
given by 

J = Zwn* (36) 


where n* is the equilibrium concentration of critical nuclei, w is the 
frequency with which a single atom joins a critical nucleus and Z' is 
a nonequilibrium (Zeldovich) factor to account for the steady-state 
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growth or decay of clusters by addition and removal of single atoms. 

The addition of single atoms to critical nucleus can occur either by 
direct adsorption from vapor or by surface diffusion of adsorbed atoms. 
However, it is believed that the surface diffusion is the more predominant 
process . 

Hirth and Pound (Reference 1-3) cite the work of Pound et al., 
which shows that surface diffusion proceeds more rapidly than direct 
adsorption from vapor by a factor of 


exp [(AF* - F*s(j)/kT] 

des 


Therefore, based on the premise that surface diffusion is the predominant 
process, co is given by the expression: 


o) = C 32 f:r*a sin G 



(37) 


where (CgpTrr* a sin 0) is the probability that an adatom is adjacent 
to a critical embryo and 



is a frequency with which an adjacent adatom will jump to join a critical 
nucleus. The nonequilibriura correction for the case of heterogeneous 
nucleation on a substrate is given by: 


z 

\3TTkTg«2/ 


(38) 


/-AF* \ 
n. = C, exp(^_^ 


(39) 


and substituting Eqs. (37), (38) and (39) in Eq. (36) and simplifying, 
we obtain: 


ma sin 0 wRT In S 
J exp 

2(o-kT/A'(e)))/2 


(-AFgd - AF*) 


kT 


c ^ 


(HO) 
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In Eq. (40), Cg is the concentration of adsorbed silicon 
concentration. Referring to Figure 2-1, silicon may be adsorbed either 
from vapor phase (step Ry) or through the decomposition of silane on 
the substrate (step Rg). Thus the concentration of adsorbed silicon 
will be determined either by the adsorption-desorption rates of silicon 
or by the heterogeneous pyrolysis of silane. 

For the adsorption of silicon ; 


Si(g) + X ►Si X 

^a “ ^a*^s^v “ ^a^s 


= kc 


CgCy 



at equilibrium, r’a = 0 and Cg = kgCgCy. If silicon is deposited through 
pyrolysis of silane, the concentration of adsorbed silicon will be 
equal to the concentration of silane molecules reacted. 


*^s = ^ra “ “ *^si “ ^h 


If the adsorption of silane, reaction of silane and desorption of hydrogen 
are all at equilibrium. 


Cgj_ = 0, = Kh2^v^h2 (molecular hydrogen adsorption) 

“ \I ^h 'J^hZ ^atomic hydrogen adsorption) 
= 0 (no hydrogen adsorption) 
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SECTION IX 
CRYSTAL GROWTH 


The mechanism of crystal growth consists of adsorption of monomer 
onto a low- index facet, surface diffusion to a ledge, association to 
the ledge, diffusion along the ledge, and incorporation into the crystal 
at a kink site. In crystal growth, ledges may be formed when two nucleated 
crystallites impinge upon one another or when a crystallite is nucleated 
in a substrate. The ideal growth rate Jq will be given by the difference 
between gross flux incident on the growing surface fp/(27rMRT) and 
the gross reevaporation flux fp /(2 ttMRT)^'^. Deviations from ideality 
are accounted for by introduction of a condensation coefficient a^. 

Thus the growth rate is given by the modified Hertz-Knudsen equation 


acf(p - p ) 

Jc = (^ 1 ) 

( 2 ttMRTs )''/2 


The growth kinetics will now be described for three different cases 
in which silicon adatoms are generated by ( 1 ) direct adsorption of 
silicon from vapor phase (step Ry), (2) dissociation of adsorbed silane 
(step Rg), or (3) dissociation of silane at ledges (step R 5 ). 


A. 


DIRECT ADSORPTION FROM VAPOR PHASE 


In this case, the continuity equation for adsorbed silicon is 
given by 


fp 


(2ttMRT.,)1/2J 



/-^des\' 

- 

N expl 1 


\ RT /J 


= D. 


8 ^Ns 


9X" 


(42) 


and the condensation coefficient is 


v/^X 


ac 


tanh 


*se 


v/2'x 


(43) 


where Ng and Ngg are the average and equilibrium concentrations of 
adsorbed silicon, A is the mean spacing between ledges, and X is the 
rms distance that an adsorbed silicon may diffuse before reevaporating, 
and given by 


X = y? 


a expi 


/^^des “ ^^sd 
2RT J 


(44) 
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At low supersaturations X is given by 


_ 19cri2 p 

X = In — r 

RTg p 


( 45 ) 


At higher supe^saturations, above a critical pressure p' given in 
Eq. (45) with X = 3^/2 X, A. is given by 


A = 3/^X 


(46) 

» 


At still higher super, saturation the nucleation of adsorbed silicon 
becomes so rapid that A'<< X and Uq becomes unity. Therefore, sub- 
stituting the proper expression for X and Eq. (43) in Eq. (42) will 
yield the expression for growth rate. 


B. DISSOCIATION OF ADSORBED SILANE 

For this case, Eq. (42) will be replaced by; 


N 


si^'s 


V exp 


^■^des\ 

."rts 7 



(47) 


Where Ng^ is the concentration of adsorbed silane and rg is the rate 
coefficient for silane decomposition. 


'^c^si’^s " 


fp 

(2TrMRTg)l/2 


(48) 


qq is still determined by Eq. (43), but (p/p') should be replaced by 
equivalent supersaturation ratio for this case. 


C. 


DISSOCIATION OF SILANE AT LEDGES 


In this case the silicon atoms produced from dissociation of 
silane are directly incorporated in the crystal, and Eq. (42) will 
be replaced by 


f pgi (2 ttMRT)-''/ 2 _ Ngi vexpi 


r-AF 


des 


92n 


RT 


s 

ax2 


(49) 
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In this equation AF(jeg, M, and D3 refer to properties of silane. 
Equations (4l) or (^3) still determine the growth rate with A referring 
to the diffusion distance for adsorbed silane. 
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SECTION X 

GLOBAL RATE OF SILICON CVD (CONCLUSIONS) 


Rate expressions for the individual steps involved in the chemical 
vapor deposition of silicon have been developed. These equations can 
be useful in two distinct areas. If the exact mechanism of silicon 
CVD and the numerical values of all the coefficients and constants 
involved in the appropriate rate equation are known, the global rate 
of silicon CVD can be determined. The rate-limiting step or steps 
can also be determined from this information. Alternatively, the rate 
expressions can bs used to evaluate the constants, determine the rate- 
limiting steps and formulate the mechanism of silicon CVD from experimental 
kinetic data on silicon CVD. The first area is modeling and simulation 
of the silicon CVD process, and the second area is basic research to 
determine the mechanism of silicon CVD. 

It is essential to know the exact mechanism of the CVD process 
to derive the global rate of silicon CVD. Published literature does 
not make it clear as to the mechanism. Any mechanism of silicon CVD 
shall consist of four basic steps: (1) mass transport of silane, 

(2) pyrolysis of silane, (3) nucleation of silicon, and (4) crystal 
growth. Rate of mass transport is uniquely given by Eq. 7- The mass 
transfer coefficient is a function of flow conditions and reactor geometry. 
Empirical correlations are available to evaluate the mass transfer 
coefficient for many different reactors. Mass transport is reported 
to be the limiting step in silicon CVD at temperatures above 1050° C. 
Pyrolysis of silane and nucleation of silicon are reported to occur 
by homogeneous and heterogeneous processes. The heterogeneous process, 
however, is believed to be the dominant process in silicon CVD. Reported 
literature indicates that surface reaction (Eq. 11) is the rate-limiting 
step in heterogeneous pyrolysis. Heterogeneous nucleation rate is 
given by Eq. (40)., Crystal growth for the heterogeneous process can 
occur either through dissociation of silane molecules adsorbed on ledges 
(Eq. 49) or through surface diffusion of adsorbed silicon atoms to 
ledges (Eq. 48). It is generally believed that growth by surface diffusion 
is the more predominant process. 

For the homogeneous process, the pyrolysis of silane is given 
by Eq. (18); the nucleation is given by Eq. (27), and crystal growth 
is given by Eq. (42). Thus we have two competing mechanisms, homogeneous 
and heterogeneoi'.G, for the chemical vapor deposition of silicon. These 
mechanisms should not be construed as real mechanisms of silicon CVD 
but only as examples to illustrate the derivation of global rate of 
silicon CVD . 

Consider the heterogeneous mechanism. The four steps involved in 
this mechanism are (1) mass transport of silane, described by rate expres- 
sion (Eq, 7), (2) heterogeneous pyrolysis with surface relation as the 
rate-limiting step (Eq. 11), (3) heterogeneous nucleation (Eq. 40), 
and (4) crystal growth by surface diffusion of adsorbed silicon (Eq. 48). 
The rates of these four steps can now be calculated for given bulk vapor 
phase concentrations of silane and hydrogen , temperatures of substrate 
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and vapor phase, and reactor pressure. Numerical values of various 
coefficients and constants are either to be taken as reported in the 
literature or estimated. Experimental data on nucleation and growth 
are lacking. If one of the steps is particularly slow (e.g., by an 
order of magnitude), that can be designated as the rate-limiting step 
and the remaining three steps can be assumed to occur at equilibrium 
conditions (i.e., highest possible rates). The rate of the limiting 
step will then be the global rate. For details of this type of analysis 
one should refer to Smith (Reference 1-2) or Hougen and Watson (Refer- 
ence 1-2), This is the same type of procedure we followed in deriving 
overall rate equations for heterogeneous pyrolysis of silane. Thus 
the two basic elements needed to use the rate expressions derived in 
this report to formulate the global rate of silicon CVD are (1) the 
exact mechanism of silicon CVD and (2) numerical values for the various 
constants involved in the rate expressions. The rate expressions will 
also be useful in analyzing experimental data on mass transport of 
silane, pyrolysis of silane and nucleation and growth of silicon crystals. 
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